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In this comment we explain the discrepancies mentioned by the authors between their results and 
ours about the influence of the gravitational quadrupole moment in the perturbative calculation of 
corrections to the precession of the periastron of quasielliptical Keplerian equatorial orbits around a 
point mass. The discrepancy appears to be consequence of two different calculations of the angular 
momentum of the orbits. 
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In [l[ the authors make use of their static and axisym- 
metric solution of Einstein's vacuum equations with a 
finite number of multipole moments [2tJ| in a system of 
coordinates adapted to multipole symmetry to derive a 
Binet equation for orbits in the equatorial plane and re- 
late it to the classical Keplerian problem. This allows 
them to write down the relativistic corrections to New- 
tonian elliptical orbits in terms of multipole moments of 
the source of the gravitational field, 
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where u is the inverse of the radial coordinate, M is 
the central mass, J is the angular momentum of the or- 
bit and V^^'^^ (u) is a generalized gravitational poten- 
tial which encloses the perturbations due to the grav- 
itational quadrupole moment Q after substracting the 
Schwarzschild and centrifugal terms. We have taken the 
gravitational constant G and the mass of the orbiting test 
particle m equal to one. 

The authors obtain a result for the angle precessed by 
the periastron in a revolution around the quasielliptical 
orbit. 
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where a is the semi-major axis of the unperturbed or- 
bital ellipse of eccentricity e and, hence, of angular mo- 
mentum J = ±-y/a(l — e'^)M and energy E — —M/2a. 
A dimensionless small parameter ( — M'^/J'^ has been 
introduced. 

In [Bl a different approach was followed to a similar 
purpose. Starting with the general metric for a stationary 
axially symmetric vacuum spacetime. 



ds^ = -f(dt - Adcj))^ + 
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where t and are the coordinates associated with the 
isometries of the spacetime and the functions /, A and 7 
depend only on the coordinates r and 9, a Binet equation 
is written for U = 1/r along timelike geodesies in the 
spacetime, 
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where E and J, respectively energy and angular momen- 
tum per unit of mass, are the conserved quantities of 
geodesic motion asociated to the isometries of the space- 
time, 

E^f{i-A4,), J = ! A{i-A^) + j^, (5) 

and the dot means derivation with respect to proper time 
along the geodesic. 

This Binet equation arises from the normalization 
condition of the velocity, v — {t,r,9,<j)), of geodesies 
parametrised by proper time, v ■ v = ~l, 
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and imposing the existence of conserved quantities in or- 
der to remove the derivatives 0, t. Binet's equation ([4]) 
is obtained dividing by 0^ and thereby eliminating the 
dependence on proper time. 

In order to compare our results with pj, we take 
Air, 0) = in order to consider only static spacetimes. 

This equation is solved perturbatively Q in powers of 
a small dimensionless parameter e = M / J and a change 
of variable — 1^0, which allows us to get rid of secular 
terms in the perturbation scheme. This frequence uj is 
responsible for the precessed angle, 

3Q ^ ,4 



UJ 



6e" 



105 



A/3 



whereas the energy is also expanded 



Ec^l + Eot^ +[-Q-10Eo^^]t\ (8) 
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and so their radius is a solution of 
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in powers of e. The term Eq is related therefore to the of the effective potential 
classical orbit. 

The discrepancy mentioned in 1] arises on identifying 
the small parameters in both expansions by means of 
C = e^. The term ZM/a in ([2]) appears to be missing in 
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The explanation is simple and is due precisely to the 
previous identification. The authors assign the classical 
Keplerian values Ec = -M/2a, = a(l - e^)M to the 
energy and momentum of the elliptical orbit. However, 
in their calculations E and J are also the conserved quan- 
tities of geodesic motion and hence J has the same mean- 
ing in both notations, though the value ~ a(l — e^)M 
should enclose the contribution of the quadrupolar mo- 
ment. 

Since we are interested in the first correction to the 
angular momentum, we perform just the classical calcu- 
lation. 

In classical mechanics the orbit of a particle moving 
under a central force of potential V{r) has two conserved 
quantities, the angular momentum J and the energy per 
unit of mass E, 



J^a + 3Q = 0. 



For > \2MQ there are two circular orbits, but we are 
interested in the exterior one. 
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since it is the one that appears as a perturbation of the 
classical orbit for small Q. 

The conserved quantities por these circular orbits, 

£; = — 5- ^, J =AfaH , 

2a , a , .a 

are seen to have simple classical corrections due to the 

presence of the quadrupole moment. 

In fact, if in ^ we include the quadrupolar correction 

to the Keplerian angular momentum, = + 3Q/a, 

the lowest order, the Schwarzschild term. 



in spherical coordinates (r, 9, 0). 

The gravitational potential due to a central mass M 
and a quadrupole M is not central, 

but considering just orbits in the equatorial plane, it acts 
as if it were central with V{r) — —M/r — Q/r^. 

For simplicity we consider circular orbits of radius 
r = a. It is clear that they cannot be used for mea- 
suring precessed angles, but still they provide an easy 
computation of ([2]) and ([7]). They are located at extrema 



6nM^ 



6t:AP 



3Q\ . 3g 



a counterterm appears that cancels out the last term in 

Hence we have shown that the apparent discrepancy 
between the formulae for the precession of the perias- 
tron of an equatorial orbit around a mass endowed with 
quadrupole moment calculated in [ij and Q is solved 
by including a first order classical contribution to the 
Keplerian angular momentum due to the gravitational 
quadrupole moment in flj. 
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